Brane-bulk energy exchange : a model with the present universe as a global attractor 
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The role of brane-bulk energy exchange and of an induced gravity term on a single braneworld of 
negative tension and vanishing effective cosmological constant is studied. It is shown that for the 
physically interesting cases of dust and radiation a unique global attractor which can realize our 
present universe (accelerating and < fl m o < 1) exists for a wide range of the parameters of the 
model. For fi m o = 0.3, independently of the other parameters, the model predicts that the equation 
of state for the dark energy today is wbe,o = —1-4, while f2 m o =0.03 leads to wbe,o = —1.03. In 
addition, during its evolution, wde crosses the m)db = -1 line to smaller values. 



I. INTRODUCTION 

In cosmologies where the present universe is realized 
as a finite point during the cosmic evolution, the answer 
to the coincidence question "why it is that today f2 m o 
and Qde,q are of the same order of magnitude" , relies on 
appropriate choice of initial conditions. By contrast, in a 
scenario in which the present universe is in its asymptotic 
era (close to a fixed point) the answer to the above ques- 
tion reduces to an appropriate choice of the parameters 
of the model. However, this latter situation is not easily 
realized if today's universe is accelerating, because: 

If the energy density of a perfect fluid with equation of 
state w>— 1/3 of any cosmological system is conserved, 
all fixed points of the system with Vl m ^ are decelerat- 
ing. 

Indeed, with p the energy density of the perfect fluid 
with conservation equation p+3(l+w)H p = 0, the Hubble 
equation of an arbitrary cosmology can be written in the 
form 



H 2 = 2 1 {p+p DE ), 



(1) 



where 7 = AttGn/S- Then, the equation governing 
Pde can always be brought into the form pus + 3(1 + 
wde)H pde = 0, where wde is time-dependent and dis- 
tiguishes one model from the other. It can be easily 
seen that d{Q, m /Q,DE)/dhia = 3(Q, m /ftDE)( w DE — w) 
and 2q = l + 'S(wfl m + WDE^DE), where fi m = 2jp/H 2 , 
Ude = ^JPde/H 2 and q= —d/aH 2 . At the fixed point 
(denoted by *) d(Q m /ilrjE)/d\na = 0. For f2 TO *^0 one 
obtains wde*~w, and 2q* = l + 3w>0. 

Thus, independently of the cosmological model, the 
only way our accelerating universe with f2 m » ^ can be 
close to a late time fixed point is by violating the standard 
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conservation equation of matter. In 4-dimensional the- 
ories, an accelerating late time cosmological phase char- 
acterized by a frozen ratio of dark matter/dark energy 
appears in coupled dark energy scenarios [l| as a result 
of the interaction of the dark matter with other energy- 
momentum components, such as scalar fields. In higher 
dimensional theories, where the universe is represented 
as a 3-brane, this violation could be the result of energy 
exchange between the brane and the bulk. In particular 
in five dimensions, a universe with fixed points character- 
ized by f2 m * 7^ 0, q* < was realized in 0] in the context 
of the Randall-Sundrum braneworld scenario with energy 
influx from the bulk. However, these fixed points cannot 
represent the present universe, since they have f2 m * > 2. 
In this paper we present a brane-bulk energy exchange 
model with induced gravity whose global attractor can 
represent today's universe. 

Let us consider an arbitrary cosmology in the form 
(|T|). Instances of such cosmologies arise in braneworld 
models or in theories with modified 4-dimensional actions 
leading to H 2 = f(p), or in cosmologies where Pde is due 
to additional fields. Assuming that as a result of some 
interaction p is not conserved, it will satisfy an equation 
of the form 



p + 3(l+w)Hp = -T. 



(2) 



Then, the equation governing pde can always be brought 
into the form 



Pde + 3{l + w DE )Hp DE = T, 



(3) 



where w EE is time and model dependent. Whenever a 
fixed point of the system satisfies 



, p = pde = 0, 



one obtains 



WDE* 



1+w 
o- 1 _ 1 ' 



(4) 



(5) 



Equation 10 is model-independent, in the sense that it 
does not depend on the form of T or the function wde (t). 
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For fi m * < 1 equation JSJ gives w db*< — 1. Specifically, 
for w = and r2 m * = ^cr>Af = 0.3 one obtains w^e* = 
— 1.4, while for r2 m * =flb ar = 0.03, wde* — — 1-03. 

The cosmology discussed in the present paper has a 
global attractor of the form @, |(5J 3j. Moreover, the 
universe during its evolution crosses the wbb = -1 bar- 
rier from higher values. This behavior is favored by 
several recent model-independent 0] as well as model- 
dependent [E IE HI analyses of the astronomical data. 



II. THE MODEL 



We consider the model described by the gravitational 
brane-bulk action 



d 5 xy^j(M 3 R-A) + / d 4 xV^h(m 2 R-V), (6) 



where R, R are the Ricci scalars of the bulk metric 
(jab and the induced metric Hab — 9ab — it-atis respec- 
tively (n A is the unit vector normal to the brane and 
A,B = 0,1,2,3,5). The bulk cosmological constant is 
A/2M 3 < 0, the brane tension is V, and the induced- 
gravity crossover scale is r c — m 2 /M 3 . 
We assume the cosmological bulk ansatz 

ds 2 = -n(t, y) 2 dt 2 + a(t, y) 2 -f l] dx l dx j + b(t, y) 2 dy 2 , (7) 

where 7y is a maximally symmetric 3-dimensional met- 
ric, parametrized by the spatial curvature k = —1,0,1. 
The non-zero components of the five-dimensional Ein- 
stein tensor are 
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where primes indicate derivatives with respect to y, while 
dots derivatives with respect to t. The five-dimensional 
Einstein equations take the usual form 

Gag = -^j^T A c\tot: (12) 

where 

Tc\tot = Tq\ v ^b + T(j\ mt B + Tc\v,b+ ?c|m,6+ Tp|i rac i(13) 

is the total energy-momentum tensor, 

diag(-A, -A, -A, -A, -A) 



Tc\v,B 



T, 



C \v,b 



diag(-V,-V,-V,-V;0) 
~- a^g{-p,p,p 7 p,0)——. 



(14) 
(15) 

(16) 



TQ\ m ,B is any possible additional energy-momentum in 
the bulk, the brane matter content T^lm^ consists of a 
perfect fluid with energy density p and pressure p, while 
the contributions arising from the scalar curvature of the 
brane are given by 



(17) 

k "')s^. (18) 



710 1 _ 
J- 1 ind — 


6m 2 


• 9 

/a 


kn 2 \ 


Kv) 


n 2 




a 2 j 


b 


Tj | ind — 


2m 2 


<a 2 


2dh 


2d 


n 2 


\a 2 


an 


a 



Assuming a Z2 symmetry around the brane, the singu- 
lar part of equations (|12f) gives the matching conditions 



P+V 
12M 3 



2n 2 



kn 2 
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(if, 



a n 



-af) (20) 



(the subscript o denotes the value on the brane), while 
from the 05, 55 components of equations l|12H we obtain 



n Q a a a b Q a Q 6M 3 



a a \a a n Q ) n 2 la a Q \a Q n a , 



a 2 



(21) 



T 55 -Ab 2 a 



6M 3 ' 
(22) 

where T05, T55 are the 05 and 55 components of TAc\m.B 
evaluated on the brane. Substituting the expressions 
(|TU|l . ipD)l in equations lj2"T)l . we obtain the semi- 

conservation law and the Raychaudhuri equation 



p + 3^(p + p) 
a„ 



2«o T o 
"67 5 



(23) 



1-- 



r 2 (p+3p-2V)] r 2 (p+3p-2V){p+V) 



-H: 



24m 2 
,2 , k 



(h 2 - 



144m 4 
r 2 ( P +V)l A- 



12m 2 



6M 3 



5 (24) 



where H Q — d /a n is the Hubble parameter of the 
brane. One can easily check that in the limit m —* 0, 
equation (|2"4"| reduces to the corresponding second order 
equation of the model without R 2\. Energy exchange 
between the brane and the bulk has also been investi- 
gated in he mm. 

Since only the 55 component of TAc\m,B enters equa- 
tion H24f> . one can derive a cosmological system that is 
largely independent of the bulk dynamics, if at the po- 
sition of the brane the contribution of this component 
relative to the bulk vacuum energy is much less impor- 
tant than the brane matter relative to the brane vacuum 
energy, or schematically 



A 



(25) 
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Then, for |A| not much larger than the Randall-Sundrum 
value V 2 /12M 3 , the term T| in equation can be 
ignored. Alternatively, the term T§ can be ignored in 
equation H24|) if simply 



A 



< 1. 



(26) 



Note that relations l|2*5|l and l|2*r?)l are only boundary con- 
ditions for Tg , which in a realistic description in terms of 
bulk fields will be translated into boundary conditions on 
these fields. In the special case where 1)25)1. I|26[) are valid 
throughout the bulk, the latter remains unperturbed by 
the exchange of energy with the brane. 

One can now check that a first integral of equation H24|) 

is 
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with x satisfying 



X + 4n H x 



rlnlT 



p+v 

6m 2 



(28) 



and T = 2T® is the discontinuity across the brane of the 
05 component of the bulk energy-momentum tensor. The 
solution of H27[) for H Q is 



2 k 



1 



6m 2 r 2 a 2 
and equation (|28ll becomes 



2(p+V) 



12 
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M3' 



-X\ > 
(29) 



X + <±n H x-- 



m 2 b n 



1± 



2\/3 



2{p+V) 12 A 



(30) 

At this point we find it convenient to employ a coor- 
dinate frame in which 6 = n = l in the above equations. 
This can be achieved by using Gauss normal coordinates 
with b(t, z) = l, and by going to the temporal gauge on 
the brane with n = l. It is also convenient to define the 
parameters 
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(31) 

(32) 

(33) 
(34) 



For a perfect fluid on the brane with equation of state 
p = wp our system is described by equations (|2*3*|) . (|2"9*|> . 



(|30|l . which simplify to (we omit the subscript o in the 
following) 



p + 3(1 + w)Hp = -T 



H 2 = fx + 2 1P ± (3^/X + 2^p+x - -5 



(35) 
(36) 



X + iH X = 24 7 r(l ± -L V /A + 24 7P + X ) , (37) 

while the second order equation (|24H for the scale factor 
becomes 

a y/X + 24 7/ o + x 



Finally, setting ?/> e 
(|38|l take the form 



VA + 24 7 p + x, equations J3HJ), 123, 



4) + 2H 



At + 2 7 p±/?V- — 
A + 6(l-3w) 7 /y 



= m- (l + 3w) 7 p±^ 



2jT 
) P 
A + 6(l-3w) 7 p 



(39) 
(40) 
(41) 



Throughout, we will assume T(p) = Ap v , with v > 0, A 
constant parameters 0, E3 • Notice that the system of 
equations l|35|) - (|37(l has the influx-outflow symmetry T — > 
—T, H — > — H, t — > — t. For T = the system reduces to 
the cosmology studied in [l4| . 

We will be referring to the upper (lower) ± solution as 
Branch A (Branch B). We shall be interested in a model 
that reduces to the Randall-Sundrum vacuum in the ab- 
sence of matter, i.e. it has vanishing effective cosmo- 
logical constant. This is achieved for /i = =f/3\/A, which, 
given that m 2 V+12M 6 is negative (positive) for branches 
A (B), is equivalent to the fine-tuning A = -F 2 /12M 3 . 
Notice that for Branch A, V is necessarily negative. Cos- 
mologies with negative brane tension in the induced grav- 
ity scenario have also been discussed in |l5j. 

Consider the case k = 0. The system possesses the 
obvious fixed point (p*, i?*, ?/>*) = (0,0, VX). However, 
for sgn(H)T < there are non-trivial fixed points, which 
are found by setting p — -0 = in equations 1)35(1 . 140JI . 
For w < 1/3 these are: 



2T(p*f 



2[i + (1 — 3u>) 7 /?* 



9(l+w) 2 p 2 

±\/9(l + w) 2 j 2 p 2 + 4/? 2 [A + 6(l-3io)7p*] (42) 
TQ>.) 



1 1 — 



(43) 



3{l+w) 

P 



7 p*V*- [A + 6(l-3io)7p»] =0. (44) 



Equation (|41|) gives 

(I). 



9(l+w) 2 p 2 ' 



(45) 
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which is positive, and also, it has the same form (as a 
function of p*) as in the absence of R. The deceleration 
parameter is found to have the value 



q* = -i, 



(46) 



which means ii* = 0. Furthermore, at this fixed point we 
find 



_2 1P * 18(1 + W ) 2 3-2, 

"m* = = ~ 7P* 



m 



A 2 



(47) 



Equation (|42|) , when expressed in terms of fl m * , has only 
one root for each branch 



j3 6(l-3w)0 ±V\(l-3w-4n^) 
27 (2n-\ + l + 3w)(n m \-l) 



(48) 



However, it can be seen from (|48f) that for — 1 < to < 1/3 
and Q m * < 1 the Branch B is inconsistent with equation 
ljI2Jl. On the contrary, Branch A with -1 < w < 1/3 
and f2 m * < 1 is consistent for < 6(1 — 3io)/3+ vA(l — 
3iu-4Q-i) < 3 v /4(l-3w) 2 /3 2 -(l + u) 2 A. Thus, since we 
are interested in realizing the present universe as a fixed 
point, Branch B should be rejected, and from now on 
we will only consider Branch A. So, we have seen until 
now that for negative brane tension, we can have a fixed 
point of our model with acceleration and < fl m * < 1. 
This behavior is qualitatively different from the one ob- 
tained in the context of the model presented in [2| (for 
— 1/3 < w < 1/3), where for positive brane tension we 
have > 2, while for negative brane tension the uni- 
verse necessarily exhibited deceleration; therefore, in that 
model the idea that the present universe is close to a fixed 
point could not be realized. 

Concerning the negative brane tension the following 
remarks are in order: (a) In the conventional, non- 
supersymmetric setting, it is well known that a negative 
tension brane with or without induced gravity is accom- 
panied by tachyonic bulk gravitational modes |l6l | ; how- 
ever, including the Gauss-Bonnet corrections relevant at 
high-energies, the tachyonic modes can be completely re- 
moved for a suitable range of the parameters 01- (b) 
As shown in 0] > m supersymmetric theories, spacetimes 
with two branes of opposite tension are stable; in particu- 
lar, there is no instability due to expanding "balooning" 
modes on the negative brane. It is, however, unclear 
what happens in models with supersymmetry unbroken 
in the bulk but softly broken on the brane. (c) Finally, 
it has been shown [jjj that with appropriate choice of 
boundary conditions, both at the linearized level as well 
as in the full theory, the gravitational potential of a mass 
on a negative tension brane has the correct 1 jr attractive 
behaviour. 



III. CRITICAL POINT ANALYSIS 

We shall restrict ourselves to the flat case fc = 0. In or- 
der to study the dynamics of the system, it is convenient 



to use (dimensionless) flatness parameters such that the 
state space is compact [2jj. Defining 



2jp_ 
D 2 



D 



D 



(49) 



where D = y/ H 2 — p, we obtain the equations 



(50) 



w m — ^ 



(i+3.)(^-i)z--^=(M^f 1 ( i-^)-^ 



-2Z(\-Z 2 ) 



2 l-Z 2 -3(l-3t«)/3V - V, 



l-Ur, 



Z' = {l-Z 2 ) (1-Z 2 ) 



l-Z 2 -3(l-3w)/? 2 /i-V, 



l-UJr, 



l + 3w 



(51) 



.(52) 



with ' = d/dT = D~ 1 d/dt. Note that -1 < Z < 1, while 
both oj's satisfy < u> < 1. The deceleration parameter 
is given by 



1- 



1 

Z 2 



l + 3w 



{1-Z' 



w m -Z 2 -3(l-3w)/3V 



1-Ulr, 



(53) 

andi?' = -HZ(q+l). The system of equations (|5l1l -ll52l 
inherits from equations (|35|I - H37(I the symmetry A —> —A, 
Z — > —Z, t — > —t. The system written in the new 
variables contains only three parameters. However, going 
back to the physical quantities H, p one will need specific 
values of two more parameters. 

It is obvious that the points with \Z\ — 1 have H = oo. 
Therefore, from (|39|l it arises that the infinite density p = 
oo big bang (big crunch) singularity, when it appears, is 
represented by one of the points with Z = 1 [Z = — l). The 
points with uj m = 1, \Z\ ^ 1, have uj' m — oo, Z' = oo and 
finite p, H; for w < 1/3, one has in addition d/a = +oo, 
i.e. divergent 4D curvature scalar on the brane. 

The system possesses, generically, the fixed point 

(a) (uj m *, Mi/,*, Z*) = (0,1,0), which corresponds to the 
fixed point (p*,^*,^*) = (0,0, V\) discussed above. 
For v < 3/2 there are in addition the fixed points 

(b) (w m *,w^,*,Z*) = (0,1,1) and (c) {uJ m *, %*) = 
(0,1,-1). All these critical points are either non- 
hyperbolic, or their characteristic matrix is not defined 
at all; thus, their stability cannot be studied by first 
order perturbation analysis. In cases like these, one 
can find non-conventional behaviors (such as saddle- 
nodes and cusps j^) of the flow-chart near the criti- 
cal points. There are two more candidate fixed points 
(d) (w m *,o;^*,^*) = (1,0,1) and (e) (w m *,w^*,^*) = 
(1,0,-1), whose existence cannot be confirmed directly 
from the dynamical system, since they make equations 
l|5ip. undetermined. Apart from the above, there 
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are other critical points given by 



2 7 



3(l+w) Z, 



(54) 



6f3 2 



(1 + 3 W KL+(1-3 W ) 1-^(1-Z 2 ) ^*-2[l-(l-Z 2 ) 2 ] 

L fl J 

= 0. (55) 

They exist only for AZ* < and correspond to the ones 
given by equations (|42() - (|44|l . For the physically interest- 
ing case w = with influx we scanned the parameter space 
and were convinced that for v^Z/2 there is always only 
one fixed point; for z^<3/2 this is an attractor (A), while 
for i/>3/2 this is a saddle (S). For w = Q, v = 3/2 there 
is either one fixed point (attractor) or no fixed points, 
depending on the parameters. For the other character- 
istic value w = 1 /3, we concluded that for v < 3/2 there 
is only one fixed point (attractor), for v>2 there is only 
one fixed point (saddle), while for 3/2 < v < 2 there are 
either two fixed points (one attractor and one saddle) 
or no fixed points at all, depending on the parameters. 
For w = 1/3, v = 3/2 there is either one fixed point (at- 
tractor) or no fixed points. Finally, for ui = 1/3, v = 2 
there is either one fixed point (saddle) or no fixed points. 
These results were obtained numerically for a wide range 
of parameters and are summarized in Tables 1 and 2. 





v < 3/2 


v = 3/2 


v > 3/2 


No. of F.P. 


1 


or 1 


1 


Nature 


A 


A 


S 



Table 1: The fixed points for w=0, influx 





i/<3/2 


i/ = 3/2 


3/2<j/<2 


v = 2 


v>2 


No. of F.P. 


1 


or 1 


or 2 


or 1 


1 


Nature 


A 


A 


A,S 


S 


S 



Table 2: The fixed points for w=l/3, influx 

The approach to an attractor described by the linear ap- 
proximation of (|51|I - H52|I is exponential in r and takes 
infinite time r for the universe to reach it. Given that 
near this fixed point the relation between the cosmic time 
t and the time r is linear, we conclude that it also takes 
infinite cosmic time to reach the attractor. 

Defining e = sgn(H), we see from (|51ll - <|52[) that the 
lines Z ~ e {y < 3/2), oj m — are orbits of the sys- 
tem. Furthermore, the family of solutions with Z e 
and dZ/duj m = Z' /uj' m w is approximately described for 
v < 3/2 by u>' m = e(l + 3w)w m (u> m — 1), and thus, they 
move away from the point (u> m *, Z*) = (1, 1), while they 
approach the point (w m *, Z„) = (1, — 1). In addition, the 
solution of this equation is u> m — [l + ce^ 14 *"'- 17 "]^ 1 , with 
c > an integration constant. Using this solution in 
equation H'/H = —Z(q+1) we find that for w = 1/3, 
H/H a = y/uj m /(l — u) m ), where H Q is another integra- 
tion constant. Then, the equation for u m (t) becomes 
dbj m /dt — — 2eu> m y / H%uj m — /Ml — u> m ) 2 , and can be inte- 
grated giving t as a function of uj m or H . Therefore, 



in the region of the big bang/big crunch singularity one 
obtains a(t) ~ yet, p(t) ~t~ 2 , as in the standard radia- 
tion dominated big-bang scenario. This means that for 
v < 3/2 the energy exchange has no observable effects 
close to the big bang/big crunch singularity. 



Z i 




FIG. 1: Influx, w — 0, u<3/2. The arrows show the direc- 
tion of increasing cosmic time. The dotted line corresponds 
to wde = — 1. The region inside (outside) the dashed line 
corresponds to acceleration (deceleration). The region with 
Z > represents expansion, while Z < represents collapse. 
The present universe is supposed to be close to the global 
attractor. 



Z i 




FIG. 2: Outflow, w = 1/3, v < 3/2. The arrows show the 
direction of increasing cosmic time. The region inside (out- 
side) the dashed line corresponds to acceleration (decelera- 
tion). The region with Z > represents expansion, while 
Z<0 represents collapse. 

Since our proposal relies on the existence of an attrac- 
tor, we shall restrict ourselves to the case ^<3/2. It is 
convenient to discuss the four possible cases separately: 
(i) w — with influx. The generic behavior of the so- 
lutions of equations (|51|I - H52|I is shown in Figure 1. We 
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see that all the expanding solutions approach the global 
attractor. Furthermore, there is a class of collapsing so- 
lutions which bounce to expanding ones. Finally, there 
are solutions which collapse all during their lifetime to 
a state with finite p and H. The physically interest- 
ing solutions are those in the upper part of the diagram 
emanating from the big bang (w, Z) w (1, 1). These solu- 
tions start with a period of deceleration. The subsequent 
evolution depends on the value of 3(3 2 /\p\, which deter- 
mines the relative position of the dashed and dotted lines. 
Specifically, for 3fl 2 /\p,\ > 1 (the case of Figure 1) one dis- 
tinguishes two possible classes of universe evolution. In 
the first, the universe crosses the dashed line entering the 
acceleration era still with wde > — 1 , and finally it crosses 
the dotted line to wde < — 1 approaching the attractor. 
In the second, while in the deceleration era, it first crosses 
the dotted line to wde < — 1, and then the dashed line 
entering the eternally accelerating era. For 30 2 /\n\ < 1, 
the dotted line lies above the dashed line, and, conse- 
quently, only the second class of trajectories exists. To 
connect with the discussion in the introduction, notice 
that the Friedmann equation i|39[l can be written in the 
form (JIJ with dark energy pde = (fiip+fi) Using (f^Ofl , 
the equation for pde takes the form Q with 



-1 



w DE - 



3(1 -W m ) 



2Z 2 -uj m ~l-6(l-3w) 



ft 2 



li 



(56) 

The global attractor ljl2"l) - (|4"4l satisfies relations Q and 
consequently, wde evolves to the value wde* given by 
(J5J. As for the bouncing solutions, they approach the 
attractor after they cross the line Z 2 — uj m , where wde 
jumps from +oo to — oo; however, the evolution of the 
observable quantities is regular. 

(ii) w = with outflow. The generic behavior in this case 
is obtained from Figure 1 by the substitution Z — > —Z 
and r — > — r, which reflects the diagram with respect to 
the uj m axis and converts attractors to repelers. 

(iii) w=l/3 with outflow. Figure 2 depicts the flow dia- 
gram of this case. Even though in the case of radiation in 
general wde > — 1/3 from equation l|56|) . there are both 
acceleration and deceleration regions. Furthermore, from 
equation JSJ it is f2 m * > 1. 

(iv) w = 1/3 with influx. This arises like in (ii) by reflec- 
tion of Figure 2 and resembles Figure 1 . 



IV. CONCLUSIONS 

In this work, we studied the role of brane-bulk en- 
ergy exchange on the cosmological evolution of a brane 



with negative tension, zero effective cosmological con- 
stant, and in the presence of the induced curvature scalar 
term in the action. Adopting the physically motivated 
p v power-law form for the energy transfer and assuming 
a cosmological constant in the bulk, an autonomous sys- 
tem of equations was isolated. In this scenario, the "dark 
energy" is a result of the geometry and the brane-bulk 
energy exchange. The negative tension of the brane is 
necessary in order to realize the present universe (accel- 
erating with 0<fi m o< 1) as being close to a future fixed 
point of the evolution equations. We studied the possi- 
ble cosmologies using bounded normalized variables and 
the corresponding global phase portraits were obtained. 
By studying the number and nature of the fixed points 
we demonstrated numericaly that our present universe 
can be easily realized as a late-time fixed point of the 
evolution. This provides an alternative answer to the co- 
incidence problem in cosmology, which does not require 
specific fine-tuning of the initial data. Furthermore, the 
equation of state for the dark energy at the attractor is 
uniquely specified by the value fl m Q. Remarkably, for 
!l m o = 0.3, one obtains Wde.o = —1-4, independently of 
the other parameters, while for the other suggestive value 
Qmo = 0.03, wde,o = —1-03. In the past, the function 
wde crosses the line wde = — 1 to larger values. 

It would be interesting to investigate if the above par- 
tial success of the present scenario persists after one tries 
to fit the supernova data and the detailed CMB spectum 
|2^ . Of course, the nature of the content of the bulk 
and of the mechanism of energy exchange with the brane 
is another crucial open question, which we hope to deal 
with in a future publication. 
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